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Abstract 

Unstable particles cannot be treated as asymptotic external states in S- 
matrix theory and when they occur as resonant states cannot be described by 
finite-order perturbation theory. The known facts concerning unstable parti- 
cles are reviewed and it is shown how to construct gauge-invariant expressions 
for matrix elements containing intermediate unstable particles and physically 
meaningful production cross-sections for unstable particles. The results and 
methodology presented are relevant for resonance physics, W^W~ and 
pair production and can be straightforwardly applied to other processes. 
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1 Introduction 



It is a fact of life that the heaviest known elementary particles, the W and bosons 
and the top quark have total decay widths that are a substantial fraction of their 
masses. Yet S'-matrix theory does not easily deal with unstable particles as they 
cannot be represented by asymptotic states. The study of width effects is timely 
and important for the theoretical understanding of current and future experimental 
results from SLC, LEP and hadron colliders. The unstable boson has by now been 
produced in vast numbers at SLC and LEP and the latter machine is soon begin pro- 
ducing W^W~ pairs. Width effects have been suggested as a possible mechanism for 
generating enhanced CP-violating observables that could provide a window into new 
physics. [0, H A number of attempts have been made to treat width effects by introduc- 
ing an imaginary part by hand into the propagator of the unstable particlep, |, ||, ||. 
Such approaches wreak havoc with unitarity and current conservation and must then 
be patched in some manner. In patching it is, in some cases, found necessary to 
introduce an imaginary part into propagators of unstable particles in the t-channel 
which is obviously incorrect. It is not easy to extend such methods to higher orders. 
A comparison of the relative merits of various approaches appears in ref.s and [Q. 
The analysis presented here suffers from none of these difficulties and is fully consis- 
tent with constraints imposed by analytic S'-matrix theory and perturbation theory. 



Much of what appears here can be found in ref.s y, |T0|, |TT1 and |112[ . 

In this paper we consider two distinct reactions involving unstable particles. The 
first is e^e~ ff near the Z^ resonance with / used to denote a generic light fermion 
species. This reaction proceeds via two distinct mechanisms. The annihilation can 
produce a Z^ boson that subsequently decays into the // pair or they can be produced 
directly without the intermediate production of a Z^ boson. For this process it is not 
too difficult to treat the reaction in terms of the stable external fermions. The Z^ 
has its 4-momentum fixed by the incoming e^e~ and therefore does not participate 
in phase space integrals. Nevertheless, to account for the Breit-Wigner resonance 
shape, some sort of resummation of the perturbation expansion is required. Done 
carelessly and without due regard to the constraints and requirements imposed by 
analytic S'-matrix theory, the resummation can lead to manifestly gauge-dependent 
results. The solution to the problem is shown to be to perform a Laurent expansion on 
the complete matrix element about the pole. This naturally generates exactly gauge- 
invariant expressions. An important point to note is that the Laurent expansion 
is much more than a mathematical trick. There is physics in the expansion. The 
leading term is the one responsible for the Breit-Wigner resonance structure and 
describes the production, propagation and subsequent decay of a physical unstable 
particle, in this case the Z". The remaining non-resonant background accounts for 
the prompt production of the final state fermions. As these processes are physically 
distinguishable in principle it is essential that their corresponding contributions to the 
matrix element be maintained separate and not combined in any way. Thus a Laurent 
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expansion should be performed even when using a perturbation expansion, such as 
the background field method |[T3[], that generates gauge-invariant Green functions. 

The second process that will be considered is the process e"'"e~ (/i/i)(/2/2) 
for energies above the Z^Z'^ production threshold. Here a complete treatment of the 
matrix element in terms of the six external stable fermions becomes unwieldy. The 
calculation would be simplified by considering just e"'"e~ —>■ Z^Z^ , that is expected to 
be the dominant source of (/i/i)(/2/2), but S'-matrix theory cannot tolerate unstable 
particles as external states. Also the invariant mass of the Z°'s is not fixed but 
varies during phase-space integrations. It will be shown, however, that from the 
matrix element for e^e~ {fifi){f2f2), it is possible to isolate a piece that describes 
propagation of Z^^s over a finite-range. This piece when suitably squared and summed 
over final state fermions yields a physically meaningful cross-section for the process 
e~^e~ Z^Z^. As such it must be exactly gauge- invariant. It turns that the required 
part of the matrix element is precisely that which is obtained by a Laurent expansion 
in the invariant masses of the Z^ bosons. 



2 Properties of Unstable Particles 



We begin by recalling what is known about unstable particles. S'-matrix theory tells 
us that unstable particles in intermediate states are associated with poles in their 
invariant momentum, s, in the S-matrix element lying off the physical sheet below 
the real s-axis. The residue of the S'-matrix element factorizes as a consequence 
of Fredholm theory (see ref. W^, p. 253). These residues can be used to define 
generalized S'-matrix elements for processes with unstable particles as external states. 
These generalized S'-matrix elements satisfy unitarity relations that are analogous to 
those for stable particles but continued off the real axisO, Unfortunately the 



unitarity relations no longer relate real quantities to one and other and it is unclear 
what they have to do with physically measurable cross-sections. 

It is known that the lifetime of an unstable particle depends on how it was pre- 
pared 



The lifetime is therefore an ill-defined concept without specifying further 
information such as that the unstable particle is in a state of definite 4-momentum. 

Veltman [0 showed that in models containing unstable particles the S'-matrix 
was unitary and causal on the Hilbert space spanned by stable particle states. That 
is to say that only stable states must be included in unitarity summations. Hence 
there is not really even room to accommodate unstable particles as external states. 
For an excellent review of some of the properties of unstable particles see ref. [ffS . 



3 Gauge Invariance near Resonance 

With the advent of LEP, the unstable Z^ has been produced in large numbers. The 
physics is described by the Standard Model Lagrangian with which calculations can 
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be done to arbitrary order in perturbation theory. The catch is that the resonance 
is a fundamentally non-perturbative object and, in order to account for the resonance 
structure, some sort of Dyson summation must be performed. As first pointed out in 
ref. the resummation can lead to a breaking of gauge-invariance. To see how this 
happens consider the process e"'"e~ — > // far from the resonance, say at PETRA 
energies of around 30 GeV. There perturbation theory works without difficulties as 
a resummation need not and should not be performed. The results generated by 
the perturbation expansion are exactly gauge-invariant in each order. Let us look 
at how the exact gauge-invariance is realized. In what follows we will use U^zzi^"^) 
to denote the transverse part of the Z^ boson one-loop self-energy. The superscript 
in parentheses indicates the loop order. The absence of a superscript will be taken 
to indicate the exact expression to all orders and used for vertices indicates the 
tree-level. The initial state vertex correction that connects the initial-state electron 
line to the Z^ is given by Viziq"^) and the final-state vertex connecting the Z° to the 
final-state fermion line by Vzfiq"^)- Although we will initially ignore the existence 
of the photon, the transverse parts of the photon self-energy and the Z-'y mixing 
are n^^(g^) and n^^(g^) and respectively. Viyiq"^) and V-yf{q^) are the initial- and 
final-state photon vertices. B{s,t) denotes one particle irreducible (IPI) corrections, 
to the matrix element. These include things like as box diagrams. 

For the process e~^e~ — > // tree-level plus one-loop corrections to the scattering 
amphtude read 

^Ml) ' (s-M2)2 ' (s - Ml) ' {s-Ml 



A{s,t) = ^ + , J,n^ + ij! + IX ' +B^^\s,t) (1) 



that is exactly gauge- invariant. Here Mz is the unphysical renormalized Z^ boson 
mass. The correction (Q) can be split into separately gauge-invariant pieces by clas- 
sifying the various contributions according to their pole structure. Writing 



nii(.) n« (M|) , ng)'(M|; 



+ nS""!^) (2) 



is - Ml) " (.-Ml) +^^^'^f^'> 
the coefficients of the pieces of (|I]) having a double, single and no pole at s = M| are 

ngk^l), (4) 

4°^ng)'(M|)yiy + V!P{MI)V^'^ + l^|V«(M|), (5) 

4°^ngr(.)i4°; + vS^%s)vi^^ + V^V^^is) + 5(^)(.,t), (6) 

respectively. Since Eq. (|1|) is exactly gauge-invariant and since the cancellation of 
gauge-dependence cannot occur between the various terms of differing pole structure 
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we must conclude that (§)-(|D are separately and exactly gauge-invariant. We will 
find that these combinations will reappear in the treatment of the resonance. 
While Eq. ([|) is gauge-invariant it blows up hopelessly as s — > M| and therefore 
cannot be used to treat the resonance. What is normally done is to perform a 
Dyson summation of the Z^ self-energy corrections to obtain an expression for the 
matrix element near resonance 

A{s, t) = ^ ' — ' + t). (7) 

s-Ml- n^l(s) 

The problem, as pointed out in ref. 0, is that A{s,t) is now gauge-dependent. 
The reason is that the Dyson summation has included all corrections of the form 

1 sf 1 each of which is gauge-dependent. In the complete matrix element, the 

gauge-dependence of these terms would be canceled by combinations of higher or- 
der self-energy corrections, n^2(s), vertex corrections, V^^^fis)-, and IPI corrections, 
B^"-\s,t). None of these are present in Eq. (0) and hence it is gauge-dependent. Far 
from resonance the gauge-dependence starts formally at 0(q;^) compared to the tree- 
level result but near the resonance where the leading s — M| in resonant denominator 
of Eq. (|^) becomes small and the gauge-dependence starts at 0{a) compared to the 
lowest order result. 

What is the solution to this problem? One that has become popular is to construct 
self-energy and vertex corrections that are gauge-invariant by themselves so that the 
resummations cannot generate spurious gauge- dependence pO|, pT] , p^] . But this is only 
symptomatic relief for a much deeper illness which the approach fails to address. In 
fact patching things in this way may generate unforeseen and undesirable side-effects 
such as shifting the pole or residue 0. 

When e"'"e~ annihilate at energies near the Z'^ resonance they produce a physical 
Z° that endures a while and then decays, normally to a fermion-antifermion pair, 
//. In a gedanken world where experimental resolution is extremely high or where 
the couplings of the Z° are extremely weak, the presence of the Z° could be detected 
as a physical particle in the way that neutrons, kaons and muons are. When e~^e~ 
annihilate they can also produce the final state // without producing a propagating 
Z^. In this case no Z° would be detected no matter how good the experimental 
resolution or how weak the couplings are made. Of course, there will always be a 
proportion of Z^ that will decay below the limit of experimental resolution. 

The // is produced by two distinguishable mechanisms and thus the matrix el- 
ement must always separate into the two corresponding pieces. How can these two 
pieces be identified in the expression for the matrix element? Consider the dressed 
propagator, for simplicity, of an unstable scalar particle. In coordinate space it is 

r d^k g-jfc-(x'-x) 
'""^^ J (27r)4 P - M2 - n(P) + te 
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The integrand has a pole at s = Sp that is a solution of the equation Sp — — Il{sp) 
0. The denominator of the integrand may be written 



(9) 



The function F{k'^) is analytic and F{sp) = (l + n'(sp)) ^. Assuming t ^ t' and 
writing Eq. (p|) in the form 



A(a;' — x) 



d^k 



-ik-{x' —x) 



integration with respect to ko yields 

d^k 



F{sp) ^F{e)-F{sp) 



k'^ - Sr. 



- Sr 



(10) 



A(x' - x) 



(27r)32A;o' 



iH^'-^)e{t' -t)+ e*^-(^'-^)^(t - t')]F{sp). (11) 



The first term in Eq. (|T^) has split into two pieces and the second has vanished because 
it has no poles. As both terms in Eq. ( [T0| ) are Lorentz invariant the vanishing of the 
second in all reference frames for t ^ t' labels it as a contact interaction. Similarly 
the first term is identifiable as an interaction of finite space-time range. It has a pole 
at k"^ = Sp and is the leading term in the Laurent expansion of the integrand about 

Sp. 

This insight can now be applied to our generic process e~^e~ ff. The exact 
scattering amplitude to all orders in perturbation theory takes the form 



A{s,t) 



V^zis)Vzfis) 



Bis,t). 



s-Ml-Ylzz{s) 

The right-hand side has a pole at s = Sp and making a Laurent gives 

Viz{sp)Fzz{sp)Vzf{sp) , Viz{s)Fzz{s)Vzf{s) - Viz{sp)Fzz{sp)Vzf{sp) 



(12) 



A{s,t) 



S - Sr, 



S - Sr 



hB{s,t). 

(13) 



where Fzz{s) is defined through the relation (||). It may be shown |]T0| that the pole 
position, residue and background are separately and exactly gauge- invariant. They 
may therefore be separately expanded as a perturbation series about M| giving 



(1)/ 



.(0) 



^(0)t.(1), 



A{s,t) 



with, to 0{a^), 



+v!^I^^is)V^ + VI^^S)V^ + vS^V^\s) + B^^\s, t) (14) 



.(0) 



rim, 



r{0) 



Sp = Ml + n« (M|) + nfzm) + nS(M|)ng'(M|). 



(15) 
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Note that the pole position, residue and background are precisely the gauge-invariant 
combinations identified in (0)-®. The overall result Eq. (|^) is clearly gauge- 
invariant and naturally separates into finite range and contact interaction as expected. 
As stated above the resonant finite range piece is associated with the production of a 
physical and the contact interaction is associated with prompt production of the 
final-state fermions via non-propagating modes of the field and box diagrams. 

It should be emphasized that nowhere in this derivation did we introduce an ad hoc 
width by hand as has been done by a number of authors ^, ||, ^. The finite width 
appeared naturally in the Laurent expansion about Sp. The subsequent expansion of 
the pole position, residue and background about the renormalized mass, Mz, is jus- 
tified because these quantities represent three independent physical observables [p!0|] . 
Since the procedure represents a well-defined sequence of expansions apphed to the 
complete matrix element, the result can be automatically guaranteed not to violate 
unitarity or gauge-invariance. 

Up to now the photon has been left out of the analysis. The Dyson summation 
will clearly be complicated by the addition of photon exchange diagrams and by Z-7 
mixing. Baulieu and Coquereaux |^ have shown how to perform the summation of 
the transverse parts exactly. The result is that the full matrix element for e"'"e~ ff 



IS 



A{s,t) = V^^is) 



n-,z(s) 



+ ^^^^'^ [s - n,,(.)] [. - Ml-Hzzis)] - ^z{s)^"'^'^ 



+ V,z{s) 



[s - U,,{s)] [s - M| - Uzz{s)] - n2 (, 



+ V,z{s) 
+ B{s,t) 



s — 11^7(5) 



[s - n,,(.)] - M| - Ylzz{s)] - lilJs 



-VzAs) 



(16) 



exactly. The first term, representing photon exchange, can be split into a resonant 
and non-resonant piece. Collecting the resonant pieces together yields 



n 



7ZI 



+ Vizis] 



n 



s-Mi- Ylzz{s) 



s — 11^7(5) 



s — 11^7(5) 



(17) 



This exact scattering amplitude has a pole at the point, Sp, satisfying the equation 



Ml - Yizzisr 



^l,{Sp) 



Sp n^^(sp) 



0. 
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Defining the function Fzz{sp) from the relation 



s-Ml- Uzzi-s) - 



Fzzis) 



[S - Sr 



(19) 



and extracting the leading term of Eq. ([T7|) in the Laurent expansion about Sp leads 
to 



A{s,t) 



Riz{sp)Rzf{sp) 



S — Sr 



^ Riz{s)Rzf{s) - Riz{Sp)Rzf{Sp) ^ Vi^{s)V^f{s) ^ ^^g) 



S - Sr. 



s — 11^7(5) 



in which 



Riz{s) = Vi^{s)- 



Tl^,z{s) 



+ V.zis) 



Fhis) 



Rzfis) 



The equation (|T7|) is exact. It is remarkably simple in structure; much more so 
than many approximate formulas that have appeared in the literature. The factor- 
ization of the residue at the pole, demanded by analytic S'-matrix theory is manifest. 
The photon exchange contribution appears in the background in a transparent form. 
As with Eq. ( |14|) the pole position, residue and background can be expanded sep- 
arately as well-behaved perturbation expansions about the renormalized mass M|. 
The background is regular in s and can be expanded as a Taylor series. The Lau- 
rent expansion has separated the finite range and contact interaction contributions to 
the matrix element and should be performed even in the case where gauge-invariant 
self-energies are produced in a consistent manner such as with the background field 
method|13|. Expanding in the way shown here, however, makes the necessity or 
advantage of gauge-invariant self-energies unclear. 

The results obtained in this section have been shown to be consistent with Ward 
identities [Ell. 



4 The Mass and Width of an Unstable Particle 

The quantity Mz in the formulas given in the previous section is the renormalized 
mass of the boson in the scheme dictated by the counterterms contained in the 
self-energies and vertex corrections. Expressions for the counterterms appearing in 
the pole position, residue and background in the Standard Model for a general renor- 
malization scheme can be found in ref. [TU|. 
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The renormalized mass is merely a bookkeeping device that has no physical con- 
tent. This is clear from the fact that, in the MS renormalization scheme, the renor- 
malized mass depends on the arbitrary scale, /x, and in the on-shell scheme it is 
gauge-dependent P3|, ^ as it is entitled to be. The physical quantity associated with 
any particle is the position of the pole of its dressed propagator. For a stable particle 
the pole position is real and it may be identified with the physical mass of the particle. 
The renormalized mass may be set equal to this physical mass by appropriate choice 
of counterterms. 

For an unstable particle the pole position is complex. It is the pole position 
as a whole that is the physically meaningful entity. The real and imaginary parts 
taken separately have no physical significance as they never appear separately in any 
expression for a physical observable. They are no more meaningful than, say, the 



modulus and argument of the pole position. Traditionally p7|, p8| for convenience the 
pole was decomposed in two possible ways, 

= M| - iTzMz (21) 

or ^ 

Sp = [Mz - ^^) . (22) 



Both definitions are arbitrary. It was pointed out independently in ref.s and |^ 



that the traditional definition of the mass lies below the on-shell renormalized mass, 
ostensibly being extracted by LEP, by 34 MeV in the case of the definition ( pij ) and 
26MeV in the case of (p2D. Sirlin||25| modified the definition of the renormalized 
mass in the on-shell scheme to be just (M^^)^ = M| — r| with Mz and T z being 
taken from (plf). It is not clear however that this can play the role of a self-consistent 
renormalized mass without violating Ward identities. Further discussion of these 



points can be found in ref. [|TT| where it is also suggested that the residue factors at 
the pole may be used to define model-independent partial widths. 

At this point a word is appropriate about the range of validity of the Laurent 
expansion used in the previous section. It is correct up to a radius determined by 
the nearest singularity to the pole, Sp. If there are thresholds in the resonance region 
then these originate branch cuts and the expansion breaks down. However some 
interesting physics then appears. In fact for the , there are a multitude thresholds 
separated, at most by the mass of a light fermion pair, lying under the resonance. 
For example, — > W^hbsc, Z° — »• 10(66) lie under the umbrella resonance but are 
sufficiently weak as make them entirely negligible. Suppose for a moment that the top 
quark had had a mass that was nit ~ Mzjl then the first-order tt threshold would lie 
under the peak of the resonance. The resonance region would be under the infiuence 
of two distinct poles, one reached by crossing the real s-axis onto the unphysical 
sheet below the threshold branch cut and the other by crossing above. Bhattacharya 
and Willenbrock|3^ have examined this scenario in a toy model and found that the 



resonance peak becomes asymmetrical being narrower on the low energy side which 
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is easily understood in that the decay width increases as the threshold is crossed. In 
the situation where mt ~ Mz/2 both poles must be associated with the Z° boson. 
The resonance then becomes a closely-spaced doublet with a separation of roughly 
100 MeV 111. 



5 Production Cross-sections for Unstable Parti- 
cles 

In the foregoing we have treated the complete process e^e~ ff. This process 
possesses a number of simplifying features. There are only a total of four exter- 
nal particles and the 4-momentum of the intermediate unstable is fixed by the 
incoming e+e~. In contrast the process e+e~ — > (/i/i)(/2/2) has an total of six ex- 
ternal particles and therefore is much more complicated to deal with especially if 
radiative corrections are to be computed. It would be advantageous if the process 
e~^e~ —>■ Z^Z^ could be treated in some meaningful way as it is expected to be the 
dominant mechanism operating in the production of the (/i/i)(/2/2)- Now instead 
of single unstable particle there are two and their 4-momenta are are no longer fixed 
but must be integrated over in phase-space integrations. The blind application of 
kinematic identities can generate complex scattering angles because of the presence 
of imaginary parts in particle self-energies. Some care is therefore required. 

As stated earlier S'-matrix theory cannot treat unstable particles as asymptotic 
states but as shown above it is possible to isolate the pieces of the matrix element for 
e+e^ (/i/i)(/2/2) that correspond to the production, propagation and subsequent 
decay of physical Z^ bosons by means of the Laurent expansion. Summing over all 
possible stable decay products then gives the total production cross-section for the 
unstable particle. 

Considering e^e" Z^Z^ avoids the need to confront the the issue of the 
Coulomb singularity that is present in e^e~ — ^ W^W~ when the W^W~ pair are 



produced at low relative velocity[pT|, |32| , [33[] . At lowest order there are no Feynman 



diagrams containing the triple vector boson vertex. 

With u and v, as usual, used to denote the spinor wavefunctions of the external 
fermions, the part of the full matrix element that can generate finite-range interactions 
is 



M =J2[Ve+TI^^Ue-]Mi{s,t,U,pl,pl) X 2 _ _ j^2 / 2^ K "^Z/i (Pl ) 



X 

P2 



where T^^, are kinematic tensors that span the tensor structure of the matrix element 
and the Mj are the associated form factors that are analytic in their arguments. It 
may be convenient, although not at all necessary, to construct the T^^, so as to be 
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individually gauge-invariant |3J]. Mj is a function of the usual Mandelstam variables 
s, t and u and of the invariant masses of the /i/i pair and pi of the /2/2- 

The expression ( p3D has been obtained by Dyson summation of the Z° self-energies. 
It is not by itself gauge-invariant as it still contains contact interaction terms that 
must be split off. In fact the complete matrix element for e^e~ (/i/i)(/2/2) 
naturally separates into four distinguishable processes according to whether not it is 
resonant or not in pi or p^. Fredholm theory guarantees the exact factorization of 
final-state contributions in the resonant channels. The four distinguishable physical 
processes are 

e+e^ — > 



e ' e 



e^e — > 

As for e~^e~ —>■ ff they become physically distinguishable in the gedanken world of 
very high detector resolution or very weak couplings. 

As it stands Eq. (^) contains the entire contribution of the matrix element for 
the first process listed above and parts of the remaining three. Additional Feynman 
diagrams must be included to compute the full matrix elements for these last three. 

Once again the Z-7 mixing has been neglected but it is a simple matter to include 
it in the manner described previously. 

Performing the Laurent expansion in pf and p^ doubly resonant part corresponding 
to the process e+e~ — >• (Z° — > /2/2) is 





> fill 








> fill 


/2/2 




/l/l 






^/2/2 


/l/l 




/2/2 





X 



Pi - Sp 

Fzzjsp) 
2 

P2 ~ "Sp 



[^hVz'f,{sp)vf,] (24) 



in which the final-state fermion currents factorize. Because Eq. (p^) is the doubly 
resonant part of the complete matrix element it must be gauge- invariant. In lowest 
it order it is 

M = h^e^T;^^e~]M,--T^[uf,V,^j,^^^^^^ (25) 
i=l Pi ~ P2 ~ 

in which T^^Mi = jf,{/pe- - /pihu/t, T^^M2 = lu{/Pe- - /P2)l^/u. The final state 
vertices then take the form V^f = ie'jp,{Pl'jL + PrIr)- As usual 7^ and •jr are the left- 
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and right-hand hehcity projection operators. The left- and right-handed couphngs of 
the Z° to a fermion / are 

_ t{ - sin^ 9wQ^ _ sin dwQ^ 

^ sin Oy/ cos Ow ' ^ cos Ow 

Squaring the matrix element and integrating over the final-state fermion momenta 
but with fixed invariant mass p\ and p\ gives the differential cross-section for physical 
production 

m\' + - + 7 + ^ - pIpI + p(pi) p(P2) 



(26) 



with 

p(p) = £E(l/?^^+l/5^^)^^^^(po)^(/) 

^ ;^-(p2-M|)2 + r|M|^^^°)^^^^ 

The ^ functions are generated in the phase space integrations and select the positive 
energy component. Thus the differential cross-section decribes the production of an 
object that has both finite range and positive energy as would be expected for a 
physical particle. The function p(p) dij)^ — Mz)9{j)q) as Q'(sp) — > that is the 
well-known result obtained by cutting in free propagator. Final state interactions 
could be built in to cross-section via the convolution kernel p{pi)-p{p2)- 

Integrating over t, pi and p| gives the total production cross-section for e+e^ —>■ 
Z^Z^ gives 

i-s /■(v^-a/pi)^ 

ct(s) = / dpi dpla{s]pl,pl) p{pl) p{pl), (27) 

JO Jo 

where 



and A = y s"^ + pf + P2 — 2spf — 2sp2 — ^pfpl- Setting pf = pi = M| reproduces the 
results of Brown and Mikaelian|]35[ . 



The result given here is superficially rather similar to those of other authors |3B, P7[ . 
The difference lies in the fact in the present case strictly only the doubly resonant 
part of the matrix element has been included and manifests itself in a change in the 
convolution kernel p. In higher orders the approach of ref.s and ||3^ will lead 
to the evaluation of off-shell matrix elements that will generally be gauge-dependent. 
In the method described here the form factors in the matrix element are evaluated 
at the pole and are therefore gauge-invariant. This also simplifies the evaluation of 
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the integrals like those in Eq. (|27|) as the higher-order form-factors are not integrated 
over. Only factors deriving from the kinematic tensors are integrated. 

After the Z°'s have propagated and decayed the final state fermions so-produced 
can interact. Because the Z°'s have propagated a finite distance from their point 
of production, the final-state interactions are expected to be suppressed. This is 
generally found to be the case[^, Such a suppression was already observed pi 



^ in the case of the e"^e~ — ^ Z — ^ // where the exact correction to the cross-section 
from the interference of initial- and final-state interference was shown to vanish near 
resonance and may be understood as suppression due the finite range propagated by 
the Z°. 

6 Conclusions 

In the foregoing paper it has been shown how to use the Laurent expansion to gen- 
erate finite-order matrix elements, for processes involving unstable particles, that are 
exactly gauge-invariant. This is achieved without the ad hoc introduction of finite 
widths and without the need for gauge-invariant self-energies or vertex corrections. 
Indeed, because of the interpretation of the resonant part as being the finite range 
interaction of an unstable particle being produced, propagating and decaying, the 
Laurent expansion should be carried out in all calculational schemes in order to pre- 
serve the separation of the matrix element into physically distinguishable parts. 

The approach was also used in obtaining physically meaningful expressions for the 
production cross-sections for the production cross-sections for unstable particles. 
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